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Methods to study event-by-event fluctuations in the 
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Theoretical calculations locate the critical point of strongly interacting matter 
(CP) at energies accessible at the CERN SPS. Event-by-event transverse momen- 
tum and multiplicity fluctuations are considered as one of the most important tools 



to search for the CP. Pilot studies of the energy dependence and the system size 
dependence of both and multiplicity fluctuations were performed by the NA49 
experiment. The NA61/SHINE ion program is a continuation of these efforts. After 

X , 

t briefly recalling the essential NA49 results on fluctuations we will discuss the tech- 

nical methods (removing Non- Target interactions) which we plan to apply for future 
Ph ■ transverse momentum and multiplicity fluctuation analyses. 



I. INTRODUCTION 

The Super Proton Synchrotron (SPS) at CERN covers one of the most interesting regions 
of the QCD phase diagram (T— //b). On the one hand there are indications (pion production, 
the kaon to pion ratio, slopes of transverse mass spectra) that the energy threshold for 
deconfinement is reached already at low SPS energies On the other hand lattice QCD 
calculations locate the critical point of strongly interacting matter in the SPS energy range 
(T CP = 162 ± 2MeV, /ig p = 360 ± 40MeV) [J. 

Fluctuations and correlations may serve as a signature of the onset of deconfinement 
(close to the phase transition the Equation of State changes rapidly which can impact 
energy dependence of fluctuations), and can can help to locate the critical point of strongly 
interacting matter (in analogy to critical opalescence we expect enlarged fluctuations close 
to the CP). For strongly interacting matter the maximum of CP signal is expected when 
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freeze-out happens close to the critical point. The position of the chemical freeze-out point 
in the (T — /xg) diagram can be varied by changing the energy and the size of the colliding 
system. 

The NA49 experiment [3] at the CERN SPS studied the energy dependence (beam ener- 
gies 20v4-158t4 GeV) and the system size dependence (p+p, C+C, Si+Si, and Pb+Pb at the 
highest SPS energy) of event-by-event transverse momentum and multiplicity fluctuations. 
There are no indications of the CP in the energy dependence of multiplicity and mean p? 
fluctuations in central Pb+Pb collisions. However, an intriguing non-monotonic dependence 
of both pt and multiplicity fluctuations was observed for the system size dependence at the 
top SPS energy [4j . Those results are consistent with a critical point located at the chemical 
freeze-out point of p+p interactions at 158^4 GeV. Quite recently a non-monotonic behavior 
of event-by-event azimuthal angle fluctuations was also observed by NA49 j^]. 



II. NA61/SHINE EXPERIMENT 

NA61/SHINE [s| is a fixed target heavy-ion experiment at the CERN SPS. Its ion pro- 
gram is a continuation of NA49 efforts. The main components of the detector, inherited from 
NA49, are four large volume time projection chambers (TPC). The Vertex TPCs (VTPC-1 
and VTPC-2), are located in the magnetic field of two super-conducting dipole magnets. 
Two other TPCs (MTPC-L and MTPC-R) are positioned downstream of the magnets sym- 
metrically to the beam line. Furthermore, the setup includes Time of Flight (TOF) walls 
and a Particle Spectator Detector (PSD). The later will be used for a very precise centrality 
determination (measurement of the energy of the projectile spectator nucleons). The excel- 
lent PSD resolution (about one nucleon in the studied energy range) is crucial especially for 
the analysis of multiplicity fluctuations. 

In the NA61/SHINE experiment hadron production in p + p, p + A, h + A, and A + A 
reactions at various energies will be analyzed. A broad experimental program is planned: 
search for the critical point, study of the properties of the onset of deconfinement, high px 
physics, and analysis of hadron spectra for the T2K neutrino experiment and for the Pierre 
Auger Observatory and KASCADE cosmic-ray experiments. Within the NA61/SHINE ion 
program we plan, for the first time in history, to perform a 2D scan with system size and 
energy. The data on p + p (2009-2010(11)), 11 B + l2 C (2011), 40 Ar + A0 Ca (2013), and 
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129 Xe + 139 La (2014) will allow to cover a broad range of the phase diagram (see Fig. [T]). 

The 2009 data taking period can be considered the beginning of the ion program in NA61. 
For each beam momentum (20, 31, 40, 80, and 158 GeV/c) we registered several million p+p 
events. Additionally, a test sample of p + p interactions at 13 GeV/c was recorded in 2010. 

III. EXTRACTION OF NON- TARGET INTERACTIONS 

For fixed target experiments such as NA61/SHINE, there is always an uncertainty whether 
the interaction occurs in the target or in some material that surrounds it. In NA61 the 
problem mostly concerns p + p interactions where the target is 20 cm long liquid hydrogen 
(LH) cylinder, and contamination may originate for example from collisions with mylar 
windows of the LH cylinder. There are numerous procedures designed to minimize the 
number of unwanted, Non- Target, interactions within data used for physics analysis. This 
includes a fine-tuned interaction trigger system as well as various quality criteria for event 
selection. With this approach one can select a group of events with the highest probability 
of coming from interactions on the target. However, for some analyses this may not be the 
most effective procedure. Event-by-event fluctuation measures can be sensitive to inclusion 
of interactions on different nuclei, in particular if this means additional p + A interactions 
within the p+p sample. Moreover, high statistics needed for those analyses can be drastically 
lowered by strict event selection criteria. 

An opposite approach is to use a more contaminated sample for analysis, but correct the 
results using data acquired from interactions on Non- Target material. Such a procedure 
is followed in NA61, where repeatedly during the data taking period the target material 
is removed. For interactions on protons the liquid hydrogen container is emptied (Target 
Empty), for interactions on solid targets, such as carbon, the target is taken out of the beam 
(Target Out). The normal experiment setup is called Target Full or Target In. 

The Target Full data sample is contaminated with Non- Target interactions. We define a 
as the fraction of Target interactions within the Full data sample: 

« = 4 (1) 
where n F is the number of Target interactions within Full data sample and n F is the total 
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number of interactions within Full data sample. For an event variable W we can write 

(W} F = -[(W) F -(l-a)(W) F NT } (2) 
a 

where (W) F is what we would like to measure in the end: event mean value for Target 
interactions within the Full data sample. (W) F is what we directly obtain from the Full 
data sample: the event mean value for the whole Full data sample. (W)h T is what we want 
to correct for: event mean value for Non- Target interactions within the Full data sample. 

In order to use Empty Target data, we assume that event mean values for Non- Target 
interactions are the same within the Full and Empty data samples: (W) F JT = {W) F T . 
Additionally, the Empty data sample consists only of Non- Target interactions (W) F T = 
{W) E , and Target interactions are only in the Full data sample (W) F = (W)t- This allows 
us to rewrite equation (j2J) using quantities that are derived directly from the acquired data 
samples: 

(W) T =-[(W) F -(l-a)(W) E \ (3) 
a 

Equation ([3]) is the main formula allowing to calculate the event mean value for interactions 
on the Target. 

Now, it is crucial to obtain the value of the a parameter. One of the possible methods 
is to analyze the distribution of the interaction vertex position v z along the beam axis. If 
one defines c as a ratio of Non- Target events within Full and Empty data samples, the a 
parameter can be expressed using c and numbers of events in Full and Empty data samples. 

F F E 

c = -> a = p (4) 

Estimation of c is possible if we assume that the distribution of v z in a given region far from 
the target is due to non target events and is identical for both Full and Empty target events 
(see Fig. [2]). The ratio of the number of events of both data samples in this region is an 
approximation of c 



n 



(5) 

v z >— 450cm 



where v z > —450 cm (the TPC gas region) describes the suggested v z range for the NA61 
experiment. 
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IV. FLUCTUATION MEASURES 



The scaled variance u is used to describe event-by-event fluctuations of the number of 
produced particles N [8] 

(X 2 ) - (AO 2 



to 



(N) 



To correct u for Non- Target interactions one needs to calculate separately 



(X) T = -[(X) F -(l-a)(X) E ) 



1 



a. 



(N 2 } T = -[(NY -(l-a)(NT\, 
a 



(6) 



(7) 



where {N) F ,(N 2 ) F , (N) E , and (N 2 ) E are quantities calculated independently for Full and 
Empty data (in the traditional way). The results are then inserted into (J6]): 



(N 2 ) T - (AQ; 

(N) T 



(8) 



The $ measure (9] is used to investigate event-by-event fluctuations of various per-particle 
quantities x, e.g. charge ($ g ), transverse momentum ($ PT ) or azimuthal angle ($,/,)• Al- 
though its formal definition contains quantities that are not event mean values as in equations 
(0 [3]), the <3> measure can also be expressed [10, lUl] through 



$ = 



(X 2 ) 2(X) (N X)_ + (X) 2 (N 2 ) 



W (xy 



\ (n) (xy 



(9) 



\ (N) (N) 2 (X) 

where X = Yli=i x i an d X2 = Z)i=i x \ are sums over particles in a given event. In this 
equation every quantity within ( ) can be calculated for Target interactions using equation 
<^fy. Two examples: 



(xy 



{(X) F -(l-a)(X) E ] 



(XX) 



T 



l(XX) F -(l-a)(XX) E ], (10) 



a a 
Using such values one can calculate $ for Target interactions, by inserting them into equation 

(ED. 



V. NA61 DATA 



The above procedure was applied to NA61 p + p data, recorded in 2009. As the data 
calibration is still in progress we used a small fraction of the available statistics and the results 
were presented at the conference only to show the method. The first glimpse of uncalibrated 
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data suggests that the correction due to contamination of Non- Target interactions will be 
small. For the scaled variance oj of the multiplicity distributions the absolute value of the 
difference between ut and u F (fluctuations obtained for Full sample only) divided by ut 
is lower than 1%. Preliminary NA61 results on the energy dependence of both transverse 
momentum and multiplicity fluctuations will be available soon. 
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FIG. 1: Two-dimensional scan of the phase diagram at SPS energies. Estimated (NA49 results 
Ref. [7]) and expected (NA61) chemical freeze-out points. 




FIG. 2: Distributions of vertex z position for p+p collisions at beam momentum 40 GeV/c. Values 
for Empty target where scaled so that the integral over v z > —450 cm is the same for Full and 
Empty setup. The enhancement around —580 cm (center of the liquid hydrogen target) corresponds 
to p+p interactions. (Uncalibrated data, small fraction of statistics, only to illustrate the method) 
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FIGURE CAPTIONS 

1. Two-dimensional scan of the phase diagram at SPS energies. Estimated (NA49 results 

h 

Ref. [7j) and expected (NA61) chemical freeze-out points. 

2. Distributions of vertex z position for p + p collisions at beam momentum 40 GeV/c. 
Values for Empty target where scaled so that the integral over v z > —450 cm is the 
same for Full and Empty setup. The enhancement around —580 cm (center of the 
liquid hydrogen target) corresponds to p + p interactions. (Uncalibrated data, small 
fraction of statistics, only to illustrate the method) 



